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ABSTRACT 

Using the results of Pichardo et al. (2005,2008), we determine regions of dynamical 
stability where planets (or discs in general) could survive in stable orbits around binary 
stellar systems. We produce this study for 161 binary stars in the Solar neighborhood 
with known orbital parameters. Additionally, we constructed numerically the discs 
(invariant loops) around five binary systems with known orbital parameters and with 
confirmed planets: HIP 10138, HIP 4954, HIP 67275, HIP 116727 and Kepler 16, as 
a test to the approximation of Pichardo ct al. (2005,2008). In each single case, the 
reported position of the planets lay within our calculated stability regions. This study 
intends to provide a guide in the search for planets around binary systems with well 
know orbital parameters, since our method defines precise limits for the stable regions, 
where discs may have established and planets formed. 

Key words: circumstellar matter, discs - binary: stars, Solar Neighborhood, exopla- 
nets. 



1 INTRODUCTION 

It is known that most low-mass main-sequence stars 
are members of binary or multiple systems (Duquennoy & 
Mayor 1991; Fisher & Marcy 1992), and in particular in the 
Solar Neighborhood, the fraction goes up to ~ 78 % (Abt 
1983). This suggests that binary formation is the primary 
branch of the star formation process (Mathieu 1994). 

Significant advances in high-angular-resolution infrared 
imaging technology have enabled large surveys of young bi- 
nary stars on a variety of star-forming regions (Mathieu et 
al. 1992, 1994). In addition, right after the discovery of the 
first extrasolar planetary system around a pulsar (Wolszc- 
zan & Frail 1992), and particularly after the first extrasolar 
planet discovered around a main sequence star (Mayor & 
Queloz 1995; Marcy & Butler 1998), observational activity 
was greatly stimulated. More recently, advances in obser- 
vational techniques and instrumentation, such as the HST 
(WFPC2 & NICMOS) imaging (Padgett et al. 1997, 1999; 
Reid et al. 2001; Borucki et al. 2010), submillimeter ima- 
ging (Smith et al. 2000), optical and infrared long-baseline 
interferometry (Quirrenback 2001a, b), millimeter and sub- 
millimeter interferometry (Launhardt et al. 2000, Launhardt 
2001, Guilloteau 2001), adaptive optics (Simon et al. 1999; 
Close 2001), spatial astrometry (Soderhjelm 1999; Quist & 
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Lindegren 2000, 2001), and microlensing (Alcock et al. 2001; 
Dong-Wook et al. 2008; Rattenbury 2009), are also availa- 
ble in binary studies. Thanks to all this new technology and 
observational work, we have now the possibility to study 
and understand better the physics of binary systems and 
the surrounding discs built during the formation stage. 

In the recent past, several planets in binary or multi- 
ple star systems have been discovered (Correia 2008; Deeg 
et al. 2008; Desidera & Barbieri 2007; Fischer et al. 2008; 
Raghavan 2006; Konacki 2005; Mugrauer et al. 2005; Eg- 
genberger et al. 2004; Sigurdsson et al. 2003; Udry et al. 
2002; Sigurdsson & Phinney 1993; Lyne 1988, etc.). Both 
suspected of formed in situ, or acquired by dynamical pro- 
cesses (Pfahl & Muterspaugh 2006). For a review of obser- 
vational techniques see Muterspaugh et al. (2010). In ad- 
dition, several recently discovered circumstellar discs whe- 
re planets are assumed to be formed, lie around close bi- 
naries (Wright et al. 2011; Prato & Weinberger 2010; De- 
sidera & Barbieri 2007; Doyle et al. 2011; Queloz et al. 
2000; Hatzes et al. 2003). In these cases the presence of 
a companion star should have a very strong influence on 
both discs and planets. From the several hundreds of ex- 
trasolar planets confirmed so far (see |http://exoplanet.eu| 
http://planetquest.jpl.nasa.gov www.exoplanets.org) to be 
around main sequence stars, about 10 % are known to reside 
in binary systems with a wide range of orbital separations. 
In almost all cases, the planet orbits in S-type configurations 
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(Dvorak 1986), while the second star acts as a perturber to 
the planetary system. A circumbinary planet (P-type or- 
bit) has recently also been detected in Kepler 16 (Doyle et 
al. 2011). This has motivated the search for stable perio- 
dic orbits around binary systems where planets (and discs 
in general) can settle down in a stable configuration. Most 
theoretical studies have focused on binaries in near-circular 
orbits (Henon 1970; Lubow & Shu 1975; Paczyriski 1977; Pa- 
paloizou & Pringle 1977; Rudak & Paczynski 1981; Bonell & 
Bastien 1992; Bate 1997; Bate & Bonnell 1997). Even very 
precise analytical methods to approximate periodic orbits in 
circular binaries are available (Nagel & Pichardo 2008). 

Due to the lack of conservation of the Jacobi integral, 
the case of eccentric binaries is qualitatively more complica- 
ted. Artymowicz & Lubow (1994) and Pichardo et al. (2005, 
2008, hereafter PSA1 and PSA2) calculate the extent of zo- 
nes in phase space available for stable, non self-intersecting 
orbits around each star and around the whole system. In 
this study we use the results of PSA1 and PSA2, to calcula- 
te stable regions for planets or discs in binary systems in the 
Solar Neighborhood with known orbital parameters such as, 
mass ratio, eccentricity and semimajor axis. 

Although the existence of stable zones, as the ones we 
are calculating here, is a necessary, but not a sufficient con- 
dition to the existence of planets (or discs in general) , if any 
stable material (planets, gas, etc.), exists in a given system, 
irrespective of their formation mechanism, they would be ne- 
cessarily located within the limits of the stable zones. In this 
direction, a fruitful line of investigation is the intersection 
between phase space available zones for the long term evolu- 
tion of planetary systems and habitable regions allowed by 
the binary system (Haghighipour et al. 2007; Haghighipour 
et al. 2010). 

We present a table with the compilation of all the bina- 
ries in the Solar Neighborhood with known orbital parame- 
ters from different sources. In the same table are presented 
the results for our calculated circumstellar and circumbinary 
stable zones around each binary system. 

This paper is organized as follows. In Section [2] we ex- 
plain briefly the method employed to calculate regions of 
stable non self-intersecting orbits around binary stars. The 
binary star sample is presented in Section [3] Results, inclu- 
ding stable regions of orbital stability for circumstellar and 
circumbinary planetary discs (and discs in general), and an 
application to observations of five binaries with observed 
planets, are given in Section 0] In Section [5] we present our 
conclusions. 



2 THE METHOD 

In the simpler case of circular binary orbits, the poten- 
tial is time-independent in the co-rotating frame and thus 
the Jacobi energy is conserved. This allows the existence of 
closed orbits that, when stable, spawn the orbital structure 
of the system (Carpintero & Aguilar 1998). In the gene- 
ral case of binaries in eccentric orbits, the problem is more 
complex, as now the potential is time-dependent. However, 
we can exploit the fact that the time-dependency is strictly 
periodic. 

A time-periodic potential in 2-D can be casted as a 3-D 
system with an autonomous hamiltonian, with the addition 



of time and the original hamiltonian as two extra dimensions 
in phase space. Regular orbits will lie on 3-D manifolds and 
be multiple periodic with three frequencies, one of which is 
given by the binary orbital frequency. If we take snapshots 
at a fixed binary phase, the projections of a regular orbit 
will lie on a 2-D manifold. If the orbit has an additional 
isolating integral of motion, this projection will now lie on 
a 1-D manifold: an invariant loop" (Maciejewski & Sparke 
1997, 2000). 

Stable invariant loops represent the generalization to 
periodically time- varying potentials of stable periodic orbits 
in steady potentials. PSA1 and PSA2, implemented a nu- 
merical method to find them. The equations of motion are 
solved in an inertial reference frame with Cartesian coordi- 
nates with the origin at the binary barycenter. An ensemble 
of test particles is launched when the binary star is at perias- 
tron, and from the line that joins both stars at that moment, 
to search for invariant loops. A more detailed explanation of 
the method and a study of the phase space in binary systems 
is in those references. In this work we employ the formulae 
from PSA1 (Eq. 6) and PSA2 (Eq. 6). These relations pro- 
vide the maximum radius of circumstellar stable zones and 
the inner radius of the circumbinary stable zone, in terms 
of the mass ratio (q = 7712/ (mi + 1712), where mi and m,2 
are the masses of the primary and secondary stars, respec- 
tively), and the eccentricity (e = \J\ — b 2 /a 2 , where a and 
6 are the semimajor and semiminor axes of the binary or- 
bit). The radius of the outer limit of the circumstellar stable 
zones from PSA1, is given by, 



Ri = Ri,E gg x [0,733(1 - e) 1 ' 20 ^ ' 07 ] 



(1) 



and a similar study in PSA2 but for the circumbinary region, 
gives a relation for the inner radius, 



Rc B (e, q) » 1,93 a (1 + l,01e ' 32 ) [q(l - q)] 



0,04:.-! 



(2) 



In these relations Ri.Egg is the approximation of Eggleton 
to the maximum radius of a circle circumscribed within the 
Roche lobe (Eggleton 1983): 



Ri,Egg/d — 



and 



0,49g; 



2/3 



1/3-, 



qi — mi/m2 = and qi = 777,2/7711 = 



1-9 



(3) 



(4) 



In Figure [T] we present a schematic figure of the geometry 
of the system. We show in this figure some of the variables 
involved in the problem. 

We must emphasize that the regions of stable orbits 
we report here, are the regions were these invariant loops 
exist, and furthermore, we restrict ourselves only to non self- 
intersecting orbits, where gas could settle and planets may 
form. It is in this sense that the term "stable region" should 
be understood in this work. Using these formulae and res- 
triction, we have calculated circumbinary and circumstellar 
radii for our sample with a total of 161 binary systems in 
the solar neighborhood with known orbital parameters. 
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Figure 1. Schematic figure of the geometry of the stable zones 
in a binary system. The black outer circle represents the inner 
boundary of the circumbinary region (eq. 01. The orange and 
green curves are the primary and secondary stars. The gray cir- 
cular areas around each star represent the circumstcllar stable 
zones (eq. [3j . The red dot at the center represents the system 
barycenter. In the upper part we show a schematic rectangle that 
we will use in Table [4] to show graphically the relative extent and 
situation of the stable regions, in relation to the stars, for each 
entry in our sample. Dotted lines join the corresponding parts in 
the upper and lower representations. 



3 THE SAMPLE 

The previous equations require, besides the stellar mass 
ratio, the semimajor axis and orbital eccentricity for each 
binary system, two parameters that are difficult to deter- 
mine observationally. There is a diversity of methods that 
have been used to estimate them. Our sample is a compi- 
lation from different sources (Jancart et al. 2005; Martin et 
al. 1998; Strigachev & Lampens 2004; Bonavita & Desidera 
2007; Holman & Wiegert 1999; Mason et al. 1999; Latham 
et al. 2002; Balera et al. 2006; Cakirli et al. 2009; Milone 
et al. 2005; Diaz et al. 2007; Konacki et al. 2010; Desidera 
& Barbieri 2007; Doyle et al. 2011). We include all binaries 
with an estimation of these parameters within 100 pc from 
the Sun (currently 161 systems). 

In Table [4] we present our sample. Columns 1 and 2, 
are the name of the systems in the Hipparchos catalogue 
and an alternate name, if it exists. Columns 3 to 6 are our 
input data: the semimajor axis, orbital eccentricity and ste- 
llar masses as reported in the references listed in column 10. 
Columns 7 to 9 list our results: the circumprimary, circums- 
tellar, and circumbinary boundaries of the stable regions. In 
the circumstellar cases, the value is the radius of the outer 
boundary. In the circumbinary, is the radius of the inner 
boundary. All lengths are given in AU and masses in so- 



lar masses. Finally, the last column is a shcematic figure 
that gives the relative sizes and positions of the stable zo- 
nes (see Figure [[J. For instance, for BD -8° 4352 (9th entry 
in the table) the circumstellar regions are symmetric and 
cover a good fraction of the inner hole of the circumbinary 
region, this is is due to the low eccentricity of the system 
and its high g(=0.5). In contrast, the binary called Gamma 
Vir (6th entry in the table), present circumstellar regions 
slightly asymmetric, due to the small mass difference bet- 
ween the stars, and quite narrow, due to its high eccentricity. 



4 RESULTS 

The presence of a stellar companion, particularly in an 
eccentric orbit, severely curtails the size and shape of the 
stable zones. While a single star has circular stable orbits 
at all radii (neglecting finite stellar size and tidal distor- 
tion effects), the presence of a stellar companion severely 
curtails the region where stable, non-self intersecting orbits 
may exist, both in extent and, to a lesser extent, shape. 

It is unclear at present the way in which these effects 
impose restrictions in the process of star formation. What 
is clear is that the effect is in the sense of inhibiting, rather 
than promote it. 

From the observational side, the statistics of the obser- 
ved systems suggests that binarity does indeed has en effect 
on planetary masses and orbits (Eggenberger et al. 2004), 
even restricting terrestrial planet formation for binary pe- 
riastron smaller than 5 AU affecting discs to within ~1 AU 
of the primary star (Quintana et al. 2007, Quintana & Lis- 
sauer 2010). 

The goal of this study is to determine the extent of 
circumstellar and circumbinary regions of stable, non-self 
intersecting orbits, as plausible regions where planets could 
have formed and may exist. It could also indicate possible 
regions of planetary formation. 

Figure [2] shows the binary semimajor axes vs. the or- 
bital eccentricity of our entire sample, split in mass ratio 
intervals. As it is already known, the region of small semi- 
major axes and high eccentricites is depleted, i.e. very close 
binaries, tidally locked in general, have eccentricities clo- 
se to zero (Duquennoy & Mayor 1991). A large percentage 
(about 60 %) of binaries in the sample have low eccentricities 
(e < 0,5), small semimajor axes (a < 50 AU) and large mass 
ratios (q > 0,4), as shown in the histograms in Figure[3] Con- 
sequently (as seen in Figure [4]), circumstellar stable regions 
in this sample have small radii (< 2 AU), with both stable 
regions (circumprimary and circumsecondary) in most sys- 
tems having similar size. On the other hand, the majority 
of circumbinary gaps have radii smaller that ~50 AU. 
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Figure 2. Binary scmimajor axes (AU) vs eccentricities of our 
sample. The mass ratio is indicated with various colors and sym- 
bols, as shown in the upper left corner. 



Figure 4. Histograms (from top to bottom): Circumprimary, cir- 
cumsecondary and circumbinary radii, for all systems in our sam- 
ple 
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Figure 3. Histograms (from top to bottom): Eccentricity, semi- 
major axis and mass ratios, for all systems in our sample. 
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4.1 Application to Real Systems 

The perturbing effect of stellar companions lead to a 
widespread belief that the presence of planets in binary sys- 
tems was very unlikely. Now we know that binaries can have 
planets, and thus should have stable regions around them 
where planet formation took place. Queloz et al. 2000, and 
Hatzes et al. 2003, discovered two giant planets in the binary 
systems GJ 86 and 7Cephei. Since then, binaries have be- 
come important targets in the search for extrasolar planets, 
particularity given their abundance. 

Until now about 70 binary systems with planets have 
been discovered (Wright et al. 2011), but only for eight of 
them, orbital parameters (semimajor axis, eccentricity and 
mass ratio) are available (Desidera & Barbieri 2007, Doyle 
et al. 2011, Muterspaugh et al. 2010). Table|4]shows our pre- 
diction for the extent of the stable orbits regions for these 
systems (in bold type letter). In all cases, where the semi- 
major axis of the planet is known, the observed planet is 
located within the predicted stable zone. We present figures 
with the calculated stable regions constructed with invariant 
loops, for the five cases where planets are confirmed, and the 
values for the semimajor axis of observed planets are known. 
In particular, notice the case of HD 120136, this is an open 
binary with a very high eccentricity, usually treated as sin- 
gle star for this reason. The large eccentricity results in a 
very narrow, circumstellar stable region. Even so, the dis- 
covered planet lies in a P-type orbit inside our predicted 
circumstellar stable region. 



4.1.1 HIP 10138 (HD 13445 or GL 86) 

This binary system is located at a distance of 10.9 ± 
0.08 pc. The companion of this object, discovered by Els et 
al. 2001, has a semimajor axis of 18.4 AU, with an eccentri- 
city of 0.4. The spectral type of the most massive compo- 
nent is Kl with a mass of 0.77 Mq and it is a white dwarf 
(Mugrauer & Neuhauser 2005). The companion has a mass 
of 0.49 Mq. At the moment, one planet was found in this 
system with a mass of M p sini = 4.0 Mj. It has a semima- 
jor axis of 0.113 AU and an eccentricity of 0.042 (Bonavita 
& Desidera 2007). For this binary the calculated stable zone 
located around each star is R ce i = 3.06 AU around the prin- 
cipal component, R ce 2 — 2.49 AU around the companion, in 
both cases this radii is the outermost radii possible to have 
stable orbits, and R c b = 58.53 AU as the innermost radii 
available for circumbinary orbits. 

In Figure [S] we present the circumprimary, circumse- 
condary and circumbinary regions of orbital stability, for 
planets in this case, to settle down, calculated at periastron 
with our method. The stellar orbits are marked in green. 
The known planet is orbiting the primary star in a very 
small orbital radius, indistinguishable in this figure. 



4.1.2 HIP 14954 (HD 19994 or 94 Cet) 

This binary system is located at a distance of 22.6 pc. 
Its semimajor axis is 120 AU and it has eccentricity of 0.26. 
The mass of the primary star is 1.35 Mq with a spectral type 
F8 V, the mass of the secondary star is 0.35 Mq. The planet 
orbiting this object has a semimajor axis of 1.428 AU with 



HD 13445 



^ 



_ dp = 0. 1 1 3 AU 
e D = 0.042 






Figure 5. Stable zones around the system HD 13445. The upper 
panel shows the circumstellar regions (primary to the left) and the 
bottom panel the circumbinary region. Notice the change in scale. 
The green lines show the stellar orbits, the barycenter is at the 
center and the system is shown at periastron. Orbital parameters: 
Mi = 0.77 M@, M 2 = 0.49 Mq, e = 0.4, a = 18.4 AU 



an eccentricity of 0.30 and M v sin i— 1.69 Mj (Desidera & 
Barbieri 2007). 

In Figure [6] circumstellar, circumbinary stable regions, 
and planetary orbit (in red) are shown. Stellar orbits (green) 
are also indicated. 



4.1.3 HIP 67275 (HD 120136 or r Boo) 

This system is located at 15.62 pc, and it has the lar- 
gest eccentricity for the cases of binary systems with known 
orbital parameters, which here are semimajor axis 245 AU, 
eccentricity of 0.91 and masses 1.35 Mq for the primary and 
0.4 Mq for the secondary component. The planet observed 
in this system has a semimajor axis of 0.048 AU while our 
approach predicts a maximum radii of 4.86 AU, the eccentri- 
city planet is 0.023 and M p sini is 4.13 (Desidera & Barbieri 
2007). 
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Figure 6. Stable zones around the system HD 19994 (94 Cet). 
The upper panel shows the circumstellar regions, primary star 
to the left secondary to the right. The bottom panel shows the 
circumbinary region, notice the change in scale. The green curves 
show the stellar orbits, the system is presented at periastron. 
Orbital parameters: M% = 1.35 M , M 2 = 0.35 M , e = 0.26, a 
= 120 AU 
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Figure 7. Stable zones around the system HD 120136 (r Boo). 
The upper panel shows the circumstellar regions, primary star 
to the left secondary to the right. The bottom panel shows the 
circumbinary region, notice the change in scale. The green curves 
show the stellar orbits, the system is presented at periastron. 
Orbital parameters: Mi = 1.35 M Q , M 2 = 0.4 M , e = 0.91, a 
= 245 AU 



In Figure we present the circumprimary, circumse- 
condary and circumbinary stable regions, for planets in this 
case, calculated at periastron. The stellar orbits are marked 
in green. The known planet is orbiting the primary star in 
a very small orbital radius, indistinguishable in this figure. 



4.1.4 HIP 116727 (HD 222404 or Gamma Cep) 

This system is located at a distance of 14.1 pc. The 
spectral type of the main component is K, and its mass is 
1.59M , while the mass of the companion is 0.4 M , the 
semimajor axis is 18.5 AU with an eccentricity of 0.36. The 
planet in this system has a semimajor axis of 2.14 AU, with 
an eccentricity 0.12 and M p sini 1.77 Mj (Desidera & Bar- 
bieri 2007). 

In Figure [8] circumstellar stable regions, circumbinary 



stable region (gray), and planetary orbit (in red) are shown. 
Stellar orbits (green) are also indicated. 

4.1.5 Kepler 16 

Recently it was found in this system a planet in cir- 
cumbinary orbit, this the first observed of this kind (on the 
circumbinary disc) . The primary star has a mass of 0.69 M , 
and the companion has a mass of 0.20 M , the semimajor 
axis is 0.22 AU with an eccentricity of 0.16 (Doyle et al. 
2011). The discovered planet has a semimajor axis of 0.71 
AU, with an eccentricity of 0.0069 and a mass of 0.33 Mj. 

In Figure [SJ circumstellar stable regions, circumbinary 
stable region, and planetary orbit (in red) are shown. Stellar 
orbits (green) are also indicated. 
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Figure 8. Stable zones around the system HD 222404 (Gamma 
Cep). The upper panel shows the circumstellar regions, primary 
star to the left secondary to the right. The bottom panel shows the 
circumbinary region, notice the change in scale. The green curves 
show the stellar orbits, the system is presented at periastron. In 
this case planet orbit is not to close to the star, line darker (red) 
around the primary star shows its orbit. Orbital parameters for 
the star: Mi = 1.59 M , M 2 = 0.4 M Q , e = 0.36, a = 18.5 AU 
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Figure 9. Stable zones around the system Kepler 16. The upper 
panel shows the circumstellar regions, primary star to the left 
secondary to the right. The bottom panel shows the circumbi- 
nary region, notice the change in scale. The green curves show 
the stellar orbits, the system is presented at periastron. Orbital 
parameters: Mi = 0.69 Mq, M 2 = 0.20 M , e = 0.16, a = 0.22 
AU 



5 CONCLUSIONS 

We have compiled a sample of binary stars with known 
orbital parameters (semimajor axes, eccentricities and ste- 
llar masses) of the Solar neighborhood and present some 
basic statistics. 

We calculate on this binary stars sample the extent of 
regions of stable non-self intersecting orbits where planets 
may exist. For this purpose, we have applied the concept 
of "invariant loops" and used the formulae of PSA1 and 
PSA2. Our approximation is ballistic, thus, the application 
is straightforward to debris discs (planets, cometary nucleii, 
asteroid belts, etc.). In the case of gas discs, further physics 



may constraint discs sizes, however, what we are providing 
here are the regions where the most important orbits of the 
binary, i.e., the ones that represent the backbone of the dy- 
namical system (the ones that are followed for the most of 
the orbits), lay. We have computed the spatial limits of these 
circumstellar and circumbinary zones for a sample of 161 bi- 
naries in the Solar neighborhood where orbital data is known 
and presented it in the form of a table where all the relevant 
parameters are provided. 

We compare our results with observations in the 5 cases 
where planets have been discovered in binary systems, and 
where semimajor axis for the planets are provided. We find 
that all the planets lay down within our computed regions of 
stability. In particular, for HD 120136, our predicted region 
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of circumstellar stability is very small, and yet the discovered 
planet lays within it. Although confrontation with a larger 
database is desirable, the current statistics is fully consistent 
with our results, proving reliable our approach. The tool 
of "invariant loops" may be very helpful in the search for 
planets in binary systems. 
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